Surface phase transitions in polydisperse hard rod fluids 
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With the simple Zwanzig model in Onsager approximation I study the effect of length poly- 
dispersity in the surface phase diagram of hard rods interacting with a hard wall. The properly 
extended interface Gibbs-Duhem equation for a polydisperse system allows us to predict the be- 
haviour of the surface tension as a function of the bulk density at the the wall-isotropic interface. 
Two groups of qualitative different bulk and surface phase diagrams are calculated from two families 
of parametrized length distribution functions p{l). This parameterization controls the law of decay 
at large Z. I also study the segregation due to polydispersity at the isotropic-nematic interface and 
the capillary nematization phenomena as a function of polydispersity. 

PACS numbers: 64.70.Md,64.75.+g,61.20.Gy 
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I. INTRODUCTION 

A molecular mixture with a finite number of compo- 
nents with different physical properties like characteris- 
tic dimension, mass or charge is a well controlled sys- 
tem in the sense that all these properties are known ex- 
actly through the study of each component in its pure 
state. Also the relative proportion of each component 
is a parameter that can be controlled in experiments. 
The situation is very different in the industrial colloidal 
systems because during its manufacture process are in- 
volved many no controlled factors at microscopic scales. 
In this case these physical properties vary in a certain 
range practically in a continuous way. Then, is more ap- 
propriate to treat these system as polydisperse in nature 
with its properties distributed in some range according to 
some continuous size distribution function. Although in 
experiments are known that many systems arc polydis- 
perse in nature there is not information about the kind of 
polydisperse distribution function that follows its phys- 
ical properties. Example of such systems are colloidal 
suspensions, polymer blends and emulsions. A common 
phenomena that have been observed experimentally in 
many polydisperse systems is that the coexisting phases 
have different particles sizes distributions and the sim- 
ulations also corroborate these results for hard spheres 

. For anisotropic particles experiments shows that the 
polydispersity makes the phase diagram much more com- 
plex including up to four different coexisting phases 
Recent experiments show that the isotropic phase of poly- 
disperse hard disk system under certain conditions can 
be more dense than the nematic phase Q . 

From theoretical point of view the polydisperse sys- 
tems composed by anisotropic particles have been in gen- 
eral modeled using a free energy functional for a binary 
mixture when the polydispersity has a clear bimodal dis- 
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tribution. The fundamental approach is then the map- 
ping of polydisperse system onto a bidisperse one. These 
bidisperse systems show similar transitions found in ex- 
periments with the organophilic polydisperse rod-like 
colloidal particles (boehmite core-particles) including a 
three phase isotropic-nematic- nematic region Never- 
theless the inclusion of bidispersity impose two relevant 
scales (the two characteristic lengths of the particles) that 
are absent in the unimodal polydisperse systems. 

Recent works include the polydispersity via a contin- 
uous size distribution function ||^. This was motivated 
by the development of the moment theory for polydis- 
perse systems which makes more accessible the theoret- 
ical study because the original infinite number of vari- 
ables that describe the whole system are reduced in this 
approximation to some moments of the distribution func- 
tion through a projection to the space extended by them 
[Q . This approximation was applied to the study of melt- 
ing in polydisperse hard spheres [^. Then, any density 
functional theory which excess part of free energy de- 
pends only on a few moments have a natural extension to 
the polydisperse case. The fundamental measure theory 
is a clear candidate for applying it to these systems. 
This functional was used in the study of the local size 
segregation of hard sphere fluid in a presence of a hard 
wall [|0|. 

Although experimentally was obtained the bulk phase 
diagram for a few polydisperse systems including also 
the solid phases there is no so many works about 
inhomogeneous polydisperse systems where the inhomo- 
geneities are created by external fields such as a hard wall 
or about the study of interfaces between the coexisting 
phases (there is a recent work ||l^ about the study of the 
surface tension and adsorption properties of the planar 
interface between coexisting fiuid phases). 

Recently have been published some works about the 
study of hard rod monodisperse system in a presence 
of a hard wall or confined in the slit showing a contin- 
uous uniaxial-biaxial nematic transition at density bel- 
low the bulk isotropic-nematic transition followed by the 
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complete wetting by nematic at the wall-isotropic inter- 
face Also was studied in these works the isotropic- 
nematic interface and the capillary nematization of the 
slit using the Zwanzig model in the Onsager approxima- 
tion to describe the hard rod system. The simulations 
and experiments fl^ confirm in the qualitative level 
all these results. 

In general, the calculation of inhomogeneous density 
profiles for mixtures of hard core anisotropic bodies us- 
ing density functional theory is a difficult problem for 
two reasons. One is the numerical problem caused by the 
presence of orientational degrees of freedom for each com- 
ponent, besides the spatial one. This complicates func- 
tional minimization. The second problem is that there 
are not so well developed and tested density functional 
for mixtures as for one-component systems. Fundamental 
measure functionals (FMF) are good candidates because 
of its natural extension to mixtures from one-component 
systems But these functionals which are very good 
at predicting the freezing of a hard spheres or the in- 
homogeneous density profiles induced by external fields 
have some difficulties if we try to extend them in 
a classical way to a system of anisotropic particles with 
continuous rotational degrees of freedom The mam 
problem here is that the Fourier transform of the direct 
correlation function calculated from the natural exten- 
sion of the FMF to anisotropic particles is isotropic at 
fc = 0, which is unphysical Nevertheless there is 

another way to construct a fundamental measure theory 
(FMT) without using the starting point of the scale parti- 
cle theory. This is the cavities method . This method 
refuses to postulate that the functional should depend 
only on certain weighted densities which are linear con- 
volutions of the density profiles. The main idea behind it 
is to ensure the dimensional cross-over of the functional 
for different configurations of cavities [|l8|. The final re- 
sult is a functional which depends on a nonlinear kernel 
that can not be decoupled in a set of one particle weighed 
densities. But this functional is impractical when there 
are rotational degrees of freedom, because the numerical 
minimization problem turns out to bee extremely hard. 

A simplification of the initial problem is to assume that 
particles are restricted to have a few discrete orienta- 
tions. Particles with different orientations are treated 
as different species and then it is completely justified to 
use the classical version of FMT for mixtures. The easi- 
est way to implement this method is the Zwanzig model 
p9[ , which describes a system of hard parallelepipeds 
with three mutual perpendicular orientations. An FMF 
for parallel hard parallelepipeds which fulfills all dimen- 
sional cross-overs 
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has been applied to the study of 
the nature of the order of nematic-smectic A transition 
I will use this functional in the Onsager limit to 
study the continuous polydisperse system. In Ref. |Q it 
has been presented the general formalism for the study 
of the bulk phase diagram of the polydisperse Zwanzig 
model in Onsager's limit. The bulk phase diagram was 
obtained for the particular choice of the Schulz size dis- 



tribution function. 

I have two purposes in the study of the polydisperse 
hard rod system. The first is that as I have already men- 
tioned the knowledge of the specific type of polydispersity 
inherent to a particular system is commonly very poor. 
Then is important to study how depends the phase dia- 
gram topology on the kind of introduced polydispersity. 
After have been obtained as general as possible a charac- 
terization of the phase behaviour for different families of 
polydispersity distribution functions (I will suppose that 
all distributions are unimodal) I will study different in- 
terfaces: the wall isotropic, the isotropic-nematic and the 
nematic polydisperse hard rod system confined between 
two walls. An interesting questions are how polydisper- 
sity affect the surface phase transitions found in and 
to study the novel surface behaviour inherent only to 
polydisperse systems. 

This paper is organized as follows. In section ^ is 
obtained formally the Onsager limit of the Fundamental 
Measure Functional for polydisperse hard parallelepipeds 
and the proper equations to calculate the bulk phase 
diagram (subsection II A ) and the different interfaces 
(subsection II B). In section ID are deduced some gen- 



eral results for thermodynamic of polydisperse interfaces, 
specifically the extension of the interface Gibbs-Duhem 
equation to polydisperse case. In section [V are pre- 
sented some important remarks about the validity of the 
thermodynamic limit in finite polydisperse systems. The 
results are presented in section]^ divided in different sub- 
sections where I study the possible bulk and surface phase 
diagra ms fo r rapidly (subsection |V A| ) and slowly (sub- 
section VB) decaying distribution functions. The last 
subsections are devoted to the study of different inter- 
faces: wall-isotropic, isotropic-nematic (subsection |V C| ) 
and a polydisperse nematic liquid crystal in the slit (sub- 
section 
section 



VI 



VP). Finally the conclusions are summarized in 



II. MODEL 

Suppose we have a multicomponent mixture of uniaxial 
parallelepipeds, all of them with the same cross section 
width CT but different lengths {L,y} {v — 1, . . . c, where c 
is the number of species). The orientation of each par- 
allelepiped is restricted to one of the three perpendic- 
ular axis X, y or z as is sketched in figure This is 
the well known Zwanzig model jl^. As I want to study 
planar interfaces I will suppose that the density profiles 
of different species are inhomogeneous only in one spa- 
tial direction and I will label them as follows: p^{L^, z), 
where /i = x,y,z labels the orientation of the species 
with length L,y. 

The excess part of free energy density according to the 
fundamental measure theory for this model is pO] 

^ex = -noHl-n-i) + + — (1) 
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FIG. 1: Zwanzig model: the uniaxial axes of parallelepipeds 
are oriented to one of the three perpendicular directions. 



where the n^'s are weighted densities, i.e. linear convo- 
lutions of the density profiles with certain weights: 



: E 

-1 }.L—x.y,z 



(2) 



where "*" means the convolution with respect to the spa- 
tial variables. These weights are characteristic functions 
whose integrals are fundamental measures of the convex 
particle |9[ |20j. The nature of the w^^^'s can be either 
scalar or vector. The scalars are 



(3) 
(4) 



where for example 51{L^) = 5{L^/2 — \z\) and 6y{L^) = 
5{a/2 — \z\). These examples allow us to recognize the 
used notation: the subindex labels the spatial coordinate 
present in the argument of the Dirac delta function. If 
the superindex and the subindex coincide then the corre- 
sponding parallelepiped edge length in the argument of 
this function is otherwise is a. The meaning of the 
labels for 9^{L^) is the same as in the preceding case, 
but 9 now represents Hevisaide function Q{x). 
The vector weighs are 



number of components, so the mixture contains paral- 
lelepipeds with lengths distributed with certain contin- 
uous function) we should replace J2u ~^ J dl in any of 
the equations that define the weighted densities, and the 
density profiles are now functions of I and z: p{l, z). We 
should remember that z is in reduced units. 

We know from Onsager that the correct variable to 
describe a nematic phase in the limit of infinite as- 
pect ratio is pL^a, so I will use in the following this 
scaled density for any specie using the same notation: 
Pfj,{l, z)L'^<j Pfj,il, z). Using the expansions of the con- 
volutions Pfj,{l, z)*S^ and Pf_i{l, z)*0i^ (/i = x, y) in powers 
of a: 



and Pi^i{l,z)a + 



(T 
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(8) 
(9) 



respectively (primes denote derivatives w.r.t. z) and ex- 
panding the excess part of the free energy density (0) in 
powers of — cf/Lq is obtained 



^$1 



(10) 



with 



$0 = ^m^^-^, (11) 
$1 - m(p + (f>o)+5I*A.TM+ n (12) 

where p — and m — (was dropped the 

dependence of and on z for notation simplicity). 
I have used the following definitions 

*M = E"^'" '^M = E^'" = (13) 



with i = 1, 2 and 



5^{L,)e^{L,)e^{L,). 



(5) 



(6) 
(7) 



where 



The component wj^xiL^) is obtained from luI^\l^) sub- 
stituting in Eq. (|) S^{L^) by 9^^{L^) and w'l?]{L^) is 
obtained from w^'^^L^) (see Eq. (^)) doing the opposite 
substitution. The rest of the components can be obtained 
using the same procedure. 

Now, suppose that we select a natural scale Lq for 
the length of the parallelepipeds, for example the mean 
length Lq = L^j^x^ (where is the molar fraction of 
the specie p). Doing the change of variable z z/Lq 
and li, = L^/Lq and taking the limit c — > oo (infinite 



p^iz) = / dlpf,{l,z), p^{z) 



m^{z) 



dUp^{l,z), m^{z) 



dlp,{l,z)^6i,{U) 

dlp:,{l,z) * 6'/, 

(15) 



with /i x, <5, ^ 6{l/2-\z\)tmAei ^Q{l/2-\z\). $o is 
then the Onsager limit whereas $1 is the first correction. 
It is interesting to notice that $0 depends only on the 
first moments of p^ {l, z) while $1 depends on the zeroth 
moments as well. 

Then, our approximation for the total free energy per 
unit area when the system is inhomogeneous in one di- 
rection is just 



jiTUajA^ Jdz ($,d( 



z) + $o(^)) 



(16) 
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where was done the proper scahng of the ideal part — > 
^idL^a which has the form 

with A|^(/) the cube of the thermal lengths of specie /i. 

A. Bulk phases 

The bulk phase diagram is calculated using the obvious 
homogeneities condition p^{l,z) — Pp(/). The total den- 
sity of the rods with length I is p{l) = PpIOj ^"^^ the 
mean density of the system p = J dlp{l). The fraction 
of rods with length I parallel to the p axis (calculated 
from P^{1) = Pfj,{l) / p{l)) is a measure of the orienta- 
tional order and so the functional $ = + <&o should 
be minimized with respect to it with the obvious con- 
straint Pfj.{l) = 1- After this functional minimization 
is obtained 



-2*„i 



(18) 



As usual, the chemical potential of rods with length / is 

calculated from (3p{l) — Using this definition and 

5p{l) 

|18|) is obtained the result 



Ppil) = In 



P{1) 



(19) 



Finally, the reduced osmotic pressure of the system is 
calculated from the thermodynamic relation: /3n = 
/ dlp{l)(3p{l) — with the following result 



/3n==p + $o. 



(20) 



The coexistence between q phases is calculated from the 
equalities between the chemical potentials and pressures 
of each type of rod for the q different phases. Select- 
ing a reference isotropic phase (referred to as the parent 
phase and labelled with superindex 0) with some fixed 
probability length distribution function then the set of 
coexisting equations are 



/i(")(0 = /i(")(0, 



n(") 



n(") 



1, 



(21) 
(22) 



The first set of equations with ( pj| ) give us the expressions 
for the densities of each phase 



(23) 



If we have an initial polydisperse parent phase with den- 
sity distribution p^^^l) then the q coexisting densities 
p'^'^^il) should obey the conservation law 



(24) 



where Xa is the volume fraction of the a phase (X^a = 
!)• 

Using (H) we obtain from (|l|)and (|2|): 

p<"Ho- l ^^b: (25) 



and using the definition p{r^(0 = P^°'\l)p'^"Hl) 

we have 



E/3 XI3J2^ 



(26) 



The definitions of the first moments, Eq. (|l^) with 
Pz{l,z) * 9i = Pz{l)l, and Eq. ( p6| ) allow us to derive 
a self consistent set of 3q equations 



TO 



(a) 



dll- 



P^x,y,z, (27) 



with a = I, . . . q. If po is the mean density of the parent 
phase pq = J dlp'^*^^{l) with p{l) its length distribution 
function, p(l) = pop{l), we have Aq unknowns: the 3q 

moments to^"\ the 9—1 independent Xa and po, and Aq 
equations to solve: the set (^) of 3g equations and the 
set of q equations for the equalities of pressures. 

Now I will specify the model for a two phase coexis- 
tence between an isotropic and a uniaxial nematic phase. 

— mi^^ = m± and 



For the nematic phase we have my 
rrix^^ = TOii^'', and for the isotropic one m^x^ 

m'P = ^TO^''^^ where m^^^ is the total mon 
isotropic phase. Defining the new variables 

2 



T = m|| + m_L — —TO 
o 



{I) 



S = TO|| — TO^, 

and using the shorthand notation 



E{1) = 3x + (1 - x)e{l), e{l) = e^^-' (e^^^ + 2) 



(28) 



(29) 



where x is the fraction of volume occupied by the 
isotropic phase, we obtain from (^7|) the equations for 
T and s: 



T = P04>T, 4>T 

The quantities p'^^\ p 



dllp{l)- 



dllpil)- 



e -2-' fe2^' - 1) 



E{1) 

-2tI ^g2si 



1 



(30) 



.(31) 



E{1) 

'^^^ (the total densities of the 



isotropic and nematic phases), and m'^^\ to*-^^ (the total 
first moments of the isotropic and nematic phases), are 
functions of r and s: 





= PO(l>plJ), 




(32) 


piN) 


— PO0p(N), 




(33) 


to(^) 


= PO<^m(J), 


(j)^(i) = 3 J dllp{l) ^^^y 


(34) 


to(^) 


= PO0m(«), 


(j>rr,(.N) = J dUp{l) ^^^y 


(35) 
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The pressures of the isotropic and nematic phases are 
related with the new variables r and s via 



/3An = Ap- 



2 (jn2 

-nv ' , 



(36) 



(r-s)(3T + s) 2mW(3T-s) 



(37) 

where AH = n(^)-n(^) and Ap = p^^'-pt-^'. To ensure 
the coexistence between the isotropic and nematic phases 
An should be equal to zero. Then, from ( ^7| ) and (p2|)- 
(p5|) is possible to calculate po as a function of t and 



(f>r)iS(t)r + 4>s) - ^4>m.(') i^'t'r " 4>s) 



(38) 



Finally, to find the coexistence between the isotropic and 
nematic phases we need to solve the closure set of two 
equations (|3^) and (^) (using expression (|3^)) with two 
unknowns r and s. These are the equations used below 
to calculate numerically, for a fixed fraction of volume 
of the isotropic phase, x, the densities of the isotropic 
and nematic coexisting phases as a function of the poly- 
dispersity. I will use the experimental nomenclature to 
classify the different coexisting curves. If there is an in- 
finitesimal amount of nematic phase coexisting with an 
isotropic phase {x = 1) the dependence of the isotropic 
density with polydispersity is called the cloud isotropic 
curve, whereas the nematic density as function of polydis- 
persity is called the shadow nematic curve. In the other 
extreme case where an infinitesimal amount of isotropic 
phase coexists with a nematic phase (x = 0) we have the 
nematic cloud and the isotropic shadow curves. 



B. Interface 

For the inhomogenous system in one spatial direc- 
tion (the z direction) was obtained (Eqs. (|l|-0) the 
Helmholtz density functional per unit area. Then the 
grand potential is, by definition. 



A 



+ ^ / / dlp^{l,z){v,,{l,z)-fi{l)), {39) 



A 



where Vfj_{l, z) is the external potential acting on a rod of 
species p and length I. A hard wall placed at z = can 
be described with the following external potential 



Vz{l,z) 




(40) 



and Vrc{l,z) = Vy{l,z) = for all z > 0. I will have 
assumed that at infinite distance from the wall the stable 



phase is the isotropic one, with mean density po, length 
distribution p{l) and mean length fixed to unity. In this 
case the chemical potential for each species is the same ( 
/i^(/) = py{l) = pzil) = p(°nO)- From ® we have 



/3p(")(0=ln(|^p(0. 



2;*(0) 



(41) 



Inserting ( ^ ) in ( p9| ) and minimizing f2 with respect to 
p^{l, z) we obtain 



Po 
3 
Po 



p{l)e 



-2iA*j,(z) 



p^x,y, 
I 



(42) 



gP(/)e-2^*^-*'^'(^)e(z-^), (43) 



^IT^ and - 



with the definition A5'p(z) *p(z 
2too/3 where ttiq = mp(oo). A self-consistent set of equa- 
tions for the first moments can be obtained from (H2), 
(^) and (llj), with the new variables ( p8|) (identifying 
TO|I with rrix and m± with my which now depend on z) 
and 



'(?(z) = 2{mz{z) - ma), 



(44) 



The result is (I drop the explicit dependence on z for 
notational simplicity) 



^po {^j dllp{l)e-^^+'^^' cosh{l 



s = -pQ I dllp{l)e-'^''+'^^^8mh{sl), 



l'Z+l/2 

dlp[l) I dz'e^ 

Jmax(z-l/2,l/2) 



In terms of these variables the equilibrium density pro- 
files, after finding the solution to (^)-(^), have the form 



Px = Y J dlpil)e' 



Py 



Po 
3 



dlp{l)e' 

2z 



(48) 
(49) 
(50) 



From ( |48| ) and (^9|) we can define an "in plane" order 
parameter as 



Qb 



Px-Py _ J dlp{l)e-'-'^+^^^ sinhisl) 
p^ + Py ~ J dlp{l)e-('»+^)^coah{sl)' 



(51) 



The model can be further simplify by a locality assump- 
tion: is approximated any convolution / * 9i{z) by f(z)l. 
Equations ( ^5| ) and (^) then remain unchanged, but ( |47| ) 
and (ISO) change to 



3P0 



2z 



dllp{l)e 



-2tI 



l<2z 

9^ = y y dlp{l)e 



-2tI 



(52) 
(53) 
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The equations (|45|)-([t6[), plus either ( [47D or (p2\), are our 
starting point for the study of the interface phase be- 
haviour of the polydisperse hard rod system in the pres- 
ence of a hard wall. 



In 1 23 it was studied, theoretically and by simulations, 
the wetting by ncmatic at the wall-isotropic interface for 
the one-component hard rod system, as well as the cap- 
illary nematization of the same system confined between 
two walls. For modeling the system the Zwanzig approxi- 
mation in the Onsager limit was used. The hard wall gen- 
erates orientational order because is entropically favored 
the alignment of rods parallel to the wall. This orienta- 
tional order enhanced by the wall translates into the den- 
sity profile in the creation of a nematic layer whose thick- 
ness increases continuously from zero to infinity when the 
bulk density is increased from Pq (the density at which 
the second order uniaxial-biaxial nematic transition oc- 
curs) to pw ( the complete wetting density) An in- 
teresting question to study is the effect of polydispersity 
on the surface phase diagram. Of course this effect de- 
pends on the length distribution function of the isotropic 
bulk phase p{l), but as I will show below, the surface 
phase behaviour for any p(l) belongs to one out of two 
qualitatively very different phase diagrams, which are de- 
termined by the behaviour of p{l) at long lengths. 

Whenever the uniaxial-biaxial nematic transition at 
the wall is a continuous transition we can calculate the 
transition density p^ for any p{l) as I will show below. 
The value of the total density and the first moment at 
contact can also be obtained analytically. 

At the wall, the value of i?(0) is — 2po/3 as can be 
easily seen from Eq. (^7|). At infinity we have d{oo) = 
because of the boundary conditions. Then we can solve 
the set of two equations ( p5| ) and ( ^ ) for any value of 
i?(z) in the interval [—2po/3,0]. 

Calling K = T + the Eqs. (^)-(^) can be put in a 
more convenient form to make these calculations: 



3P0 



dllp{l)e~'^^ smhsl, 



3P0 



dllp{l) 



coshs^ — 1) 



(54) 
(55) 



The variable s calculated from the equation s — f{K,,s), 
where / is the right hand side of Eq. (jsj), is a kind of 
surface order parameter, which is zero below the tran- 
sition density Pq and is positive above it (I will assume 
that the symmetry breaking at the wall occurs in the x 
direction). 

If the transition is continuous then f'^^K* , 0) = 1 should 
be obey at the transition point. Using this fact and com- 
bining Eqs. (^) and (|55|), we obtain one single equation 
to calculate the bifurcation value k*. 



dllp{l)e~'''^ (1 - K*l) ^ 2, 



and from it the bifurcation density. 



dlfp{l)e-'' 



(56) 



(57) 



Equations (|5^) and ( pTj ) will be used below for the calcu- 
lations of the second order uniaxial-biaxial nematic tran- 
sition at the wall for different length distribution func- 
tions. 



III. THERMODYNAMICS OF POLYDISPERSE 
INTERFACES 

Once the grand potential is minimized we can calculate 
the surface tension 



7 ^ 



(58) 



where L is the length of interface and Hq the bulk os- 
motic pressure. To characterize surface transitions (as 
wetting) is important to calculate the density adsorption 
coefficient. For the polydisperse case I will show that, 
apart from it, the other moment adsorption coefficients 
of the total density distribution function are also rele- 
vant. In the excess part of free energy functional of this 
model only the first moment appears ,so the quantities 
we are referring to are 



Pa 



(59) 



r^'^Po) = V / dz(m^(2)-mo). (60) 
,, Ja 



In a mixture of a finite number of species the interface 
Gibbs-Duhem equation relates surface tension changes 
with changes of the chemical potentials of the different 
components as in 



(61) 



where F^ is the adsorption coefficient of component i^, 
and pi, its chemical potential. The natural extension of 
Eq. (|6l| ) to the polydisperse case is 



dl = - j dlT{l)dp{l), 



(62) 



where r{l) = dz {p{l , z) - p<-°\l , 00)) . Supp ose the 
system is composed by a continuous mixture of asymmet- 
ric particles with one of its characteristic lengths, say I 
(from the complete set of characteristic lengths that de- 
fine the geometry of the particle), taking different val- 
ues according to some fixed length distribution function 
p(Z), whereas the other lengths are constant. In general, 
the chemical potential of the species with length in 
the bulk isotropic phase with total density po can be ex- 
pressed, according to the scale particle theory, as 



D' 



p{l) = In (poP(O) + «fe(/'o)^* 



(63) 



fc=0 
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where the first term comes from the ideal part and the 
second one from the excess part of the free energy density, 
which in any scale particle theory is a polynomial on I 
of degree D' with coefficients depending on po- D' is a 
number between 1 and D (the dimension of the system) 
which depends on the specific geometry of particles. For 
example, in three dimensions we have D' — S for hard 
spheres, D' = 2 for uniaxial oblate parallelepipeds, or 
D' = 1 for uniaxial prolate parallelepipeds. 

A change in the bulk density po implies a change in 
the surface tension, and using Eqs. (^) and ( |63| ) we can 
find an expression for the derivative of the surface tension 
with respect to the bulk density. 



dpo 



r(o) 

Po 



D' 



k=0 



where 



r 



(fe) 



(64) 



(65) 



is the fc- moment adsorption coefficient. The fc- moment 
is defined as usual: M'-'^^z) = J dll''p{l,z). 

This important result tells us that the slope of the sur- 
face tension depends in general on different moment ad- 
sorption coefficients of the density distribution function 
p{l, z). It is thus possible to find a qualitatively different 
behaviour of the surface tension as compared to that of 
the one-component system, for which the slope depends 
only on the density adsorption coefficient, as we can see 
after substituting p{l, z) = p{z)6{lo — I) in (|6^): 



c?7 
dpo 



Po 



k=0 



r(o). 



(66) 



For example, as I will show below, fixing the degree 
of polydispersity we can find some range of po's for 
which the density adsorption coefficient r^"-* and the 
slope d'y/dpo are both negative, in contrast with the 
mono-disperse case behaviour (see Eq. (^6|)), for which 
d'y/dpo and F^"-* have always opposite signs. This is due 
to the fact that the first moment adsorption coefficient 
F^-'^) is still positive as a consequence of the higher deple- 
tion effect that the wall induces on the large rods. Then, 
from Eq. ( |6^ ) we can see that the final sign depends 
on the modules of F(°)/po and a[{po)r('^^ (D' ^ 1 and 
o-o{po) — for this model). 

For the sake of further illustrating the validity of Eq. 
(U) I will deduce the surface tension as a function of 
density po near the complete wetting transition, po ~ Pw 
The adsorption coefficients diverge logarithmically near 
the transition point, according to 



fW ^i^fclni + CTfc, 



(67) 



where t = 1 — po/pw- Considering the Zwanzig model for 
prolate parallelepipeds in the Onsager approximation, for 



which all are zero except ai = 4po/3, and inserting 
(|^) into (^) we have 

(— + -1^1) hit -cTi. (68) 

dpo Po o Po o 

Integrating (|6^) with respect to po we obtain 

A7 = iyQ£{t) - ao ln(l -t)+t [£>i(ln i - 1) + cti] , (69) 

where A7 — "f — jw, with 7^ the surface tension value at 
the transition point, j>i = 4/5^^i/3, ai = 4p„(Ti/3, and 
C{t) = J^dx{l-x)''^\nx. 

I will use Eq. ( |69| ) for testing the validity of (|^ with 
the exact surface tension calculations near the wetting 
transition. 



IV. POLYDISPERSE SYSTEMS WITH FINITE 
NUMBER OF PARTICLES 

In real systems the total number of particles although 
large is a finite number. For one-component systems is 
completely justify to take the thermodynamic limit but 
this is not the case for polydisperse systems. The funda- 
mental reason of this difference is that the finite systems 
have only a few species with characteristic lengths dis- 
tributed in the tail (far from the mean value) of the size 
distribution function. The majority of theoretical calcu- 
lations on polydisperse systems have in common the hy- 
pothesis that for any length there are an infinite number 
of particles (in the thermodynamic sense) distributed ac- 
cording to the so called parent size distribution function 
p{l) (normally infinite ranged). Although this hypothe- 
sis is in general not justified the majority of the obtained 
results reflect qualitatively the experimental behaviour 
because normally the decaying law of p{l) for large Vs 
is so rapidly that as is shown here the particles with 
large lengths don't affect too much the properties of the 
different coexisting phases. For low decaying p(Z)'s like 
the log-normal distribution the species with large lengths 
have a strong influence in the phase behaviour of the sys- 
tem. Then, the phase diagrams obtained through theo- 
retical calculations using the preceding hypothesis can 
differ very much from the real ones. For this reason be- 
comes crucial to estimate the number of particles present 
in the tail of p{l) and if is highly enough to justify the 
thermodynamic limit. 

Suppose we have a polydisperse inflnite system de- 
scribed with infinite ranged size distribution function p(^) 
and that we take a sample of N particles (with very 
large) which coincide with our real system. Then, there 
is a maximum extreme value Lmax for the lengths of 
particles which obviously depend on N and on the de- 
gree of polydispersity. I will obtain this relation using 
the asymptotic properties of sampling distributions (for 
more information see the reference [p4|). Arranging the 
A^ sample values in order of magnitude and considering 
the A't'th value from the top, for Nt = 1 we obtain the 
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extreme value. If I denote the length of the iVt 'th specie 
from the top, the probability element that among the 
N sample values, N — Nt are less than I, and Nt — 1 
are greater than I + dl, while the remaining value falls 
between I and I + dl is: 



gj,,{l)dl ^ NCZ,Z\F{l f-''^ {1 



F{l)f-'p{l)dl, {70) 



where F(l) = JjJ dtp{t) and C]^Zi is the combinatorial 
number. Introducing the new variable ^ = N{1 — F{1)) 
we have that the distribution function that follow this 
variable in the asymptotic limit of large A^'s is hN^{^) = 
^Wt-i gxp (_^)/p(7Vf ), where r(x) is the Gamma func- 
tion. The variable / depends on ^ and N by the relation 



(71) 



^:^N p{t)dt. 



If the asymptotic solution of this equation is Z = /jv(0 + 
0{u{N)) (where u{N) is the next term in the asymptotic 
expansion) we can estimate the mean length (the average 
is taken over an infinite number of samples with fixed N) 
of the A^t'th specie from the top as 



LNiNt,A) 



d^hN,iOfNiO+Oiu{N)), (72) 



where was shown explicitly the dependence oi on 
the degree of polydispersity A. Fixing the total number 
N , the maximun extreme value is obtained as Lmax — 
LAr(l, A). Then, the number of particles with lengths 
in the interval [(1 — r)L,„ax, imax] (with r < 1) is the 
solution of the equation 



ijv(^t,A), 



(73) 



with respect to Nt. Here I will use two different fam- 
ilies of length distribution functions. All them have 
in common that the mean rod length is fixed to one 
(/ dllp{l) = 1). Is imposed also that hm;_op(0 = 0; 
i.e. the probability of finding rods with zero length is 
zero. As a measure of the degree of polydispersity I se- 
lect the parameter A = \/a^ — 1, where a is the mean 
deviation = J dlPp{l). The last constraint that I im- 
pose to p{l) is limA->oP(0 — ^(.^ ~ 1); so is recovered the 
one-component system in the limit of small polydispersi- 
ties. 

The first one is defined as 

^(^^ = r[(J+i)Vi] ^'^^"P["^^--^^^'J' ^^^^ 
_ r[(^ + 2)g-i] 
- r[(. + i)g-r 

The number q > I controls the decaying law at large ^'s, 
and v > controls the degree of polydispersity which 
maximun value 



Aniax 



/ r(g-i)r(3g- 

r2(2g-i) 



(75) 



is reached at = 0. When g = 1 (74) coincide with the 
Schulz's distribution function for which Amax = 1- 
The second family is 



p(0 = C„,,expha|ln(Z//o)n 



(76) 



where Ca^q and Iq depend on a and q and are calculated 
from the conditions / dlp{l) = J dllp{l) — 1. Specially 
for q = 2 we have the log-normal distribution: 



P{1) 



In' 



^27rln(l + A2) 



exp 



(1 + A2)l; 



21n(l + A2) 



,(77) 



where < A < cxd. Applying the described formalism to 
the distribution function family (^3) we obtain 



LNiNuA) 



1/9 



1 - rfo ( 1 + y + (g - l)rfi ) , (78) 



where were neglected terms of order ^ In^^'' "^{N) and I 
have defined =^ (ly + l)/q - I, t ^ ln[7V/r(/i + 1)] and 



dk = — 
qt 



1 flk{Nt) 



2k 



(79) 



with 7a:(x) the Polygamma functions of order k. 
For the log-normal p(Z) (see Eq. |7^) we obtain 

(47rt)'/('^^r(iVt) 



where were neglected terms of order of unity. Now 
t = IniV and A = y^2t ln"^(l + A2). Fixing N = lO^^ 

and r = 4/7 I solved the equation ([73| ) for Nt using Ln 
from (^) with g = 1,2 and from (g^. The results {Nt 
as a function of A) are shown in figure ^. As we can 
see from the results for low decaying distributions the 
tails contain less number of particles. This behaviour can 
be understood if we take into account that fixing N the 
probability to find large particles strongly decrease with 
q and then Lmax also decrease in so way that the number 
of particles present in the tail (with length rimax) is an 
increasing function of q. Another interesting feature is 
that Nt decrease with A (except for g = 1 and A < 0.3) 
which can be explained by the same reason because Lmax 
is an increasing function of A. The main conclusion that 
we can extract from these results are that the thermo- 
dynamic limit for the species belonging to the tail of the 
length distribution function is less justified from low de- 
caying p{iys and for large polydispersities. 

Following the usual theoretical works on polydisperse 
systems I will fix a parent p{l) but taking into account 
that the system contains a finite number of species will be 
used a truncated p{l). The main purpose is to estimate 
the number of particles present in the tails. Following the 
already described procedure I fix now Lmax and solving 



equation Ljv(l,A) 



for N I find the dependence 
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A 



FIG. 2: Number of species with lengths in the interval 
] (r = 4/7) as a function of polydispersity 
for fixed total number of particles: = 10^^. With solid, 
dashed and dotted lines are shown the results for log-normal, 
Schultz's distribution function and p{l) from (j?^) with q = 2 
respectively. The point represent the value Amax for q — 2. 
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A 

FIG. 3: Number of species with lengths in the interval [(1 — 
r)Lmax,imax] {r = 4/7) as a function of polydispersity for 
fixed Lniax = 20, 35 and 50 (the dotted, solid and dashed 
lines respectively). 

iV(A). Obviously each point in the p — A phase diagram 
now represent systems with different number of parti- 
cles. From equation ( [73| ) I calculate A^f(A) for the same 
r — A/7 which is shown in figure ^ for the log-normal 
distribution with imax = 20, 35 and 50. Although the 
relative fraction Nt/N is an increasing function of A 
as should be, Nt decrease with polydispersity. Taking 
into account that the shadow phase fill an infinitesimal 
amount of the total volume, the number of its particles 
is less in orders of magnitude compared with the cloud 
phase. Then the fractionation of Nt ^ 10'^ particles from 
the parent to the shadow phase can change dramatically 
the mean length. We can see from figure || that until 
A ~ 0.6 we have Nt ^ 10^ so the results presented in 
the next sections for the log-normal distribution function 
should be qualitatively correct in this range of polydis- 
persity values. 



V. RESULTS 

The important quantities to characterize the surface 
behaviour of our hard rod system interacting with hard 
wall are the fc-moment adsorption coefficients {k = 0, 1) 



and the surface tension (see Eqs. (p8|), ( p9| ) and ( |60|)) . 
This quantities are calculated after minimization of (|39|) 
to obtain the equilibrium density profiles Pfj,{z) (see Eqs. 
(^8|)- (^o|)). However, we can gain insight into the surface 
phase diagram for this model without carrying out the 
full minimization as I will show bellow. 

Because the hard wall does not allow a rod with length 
I perpendicular to it to be at a distance less than 1/2, we 
have mz(0) = 0, or equivalently 'd{0) = —2po/3. Then 
Eqs. ( ^ and ( ^ ) can be solved to obtain t(0) and s(0), 
which yield inserted in (|8|) and ( |49| ) PxiO) and Py{0). 
Finally, p"^ = J dlpz{lJ/2), the total contact density of 
rods perpendicular to the wall, is calculated using the 
contact value theorem valid for any mixture of particles 
interacting via hard core. This theorem for a mixture of 
hard core particles reads 

f3Il = Y.PM^J)- (81) 

The extension of this theorem to the polydisperse 
Zwanzig model in the Onsager limit is 

E f dlp,{l,l,)^pAO) + Py{0) + P'z, {82) 

^—x^y.z 

where /3n = po(l + 2po/3) is the bulk pressure and Ix — 
ly = 0, Iz = 1/2. Finally, it is necessary to know for 
which value of the bulk isotropic density po there appears 
a macroscopically thick layer of the nematic phase at the 
wall-isotropic interface. For complete wetting this value 
can be calculated from ( ^0| ) and (|3^) with x = 1 (the 
fraction of volume occupied by the coexisting nematic 
phase is zero). 

Having all this information about the surface structure 
we can determine the possible scenarios of surface phase 
transitions occurring before the wetting transition, de- 
pending on the form of the length distribution function 
p{l). For example, we can distinguish the nature of the 
surface transition (second or first order). For the former 
we obtain the exact value of the transition point us- 
ing equation (p7|), and for the latter we can also calculate 
the exact value if we use the local model, and produce an 
estimation for the nonlocal model, as we will show later. 

I have found that the surface phase behaviour depends 
on the decaying law at large lengths of the length dis- 
tribution function p(l). I have found two qualitatively 
different types of surface phase diagrams, one of them 
coming from p(/)'s decaying exponentially or faster, and 
the other from slowly decaying p(/)'s. 



A. Bulk and surface phase diagrams for rapidly 
decaying distribution functions 

I begin the study with p{l) from the first family 
(fz^). Taking into account that I pretend to study sys- 
tems with very large but finite number of particles (see 
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the discussion done in the preceding section) the cor- 
rect way to do this is to select the length distribution 
function equal to p(Z)6(Li„ax — I), where Lmax is the 
maximun extreme value for the lengths of particles and 
0(x) is the Hevisaide function. With the constraints 
jj'""'" dll''p{l) = 1 (fc = 0, 1) the equation (|7|) keeps the 
same but substituting r(a;) by the Incomplete Gamma 
function r[(c,y_gL,nax)'', a;] where c^^q is now the solution 
of the equation 



r[Kgimax)«,(!^+l)g-l] 



(83) 



Using Eqs. ( |56| ) and I have calculated the second 
order uniaxial-biaxial nematic transition (or isotropic- 
nematic transition in the plane of the wall), i.e. the bulk 
density value pb at which the solution nixifi) ^ '7iy(0) 
bifurcates, using this family of p(/)'s for several values of 
q. When the bulk density po is greater than this value we 
find a phase with preferential alignment of rods parallel 
to the X axis, i.e. > Py,Pz with py ^ p^ due to the 
presence of the wall. This is the reason why in this 
phase was called biaxial nematic phase. All results us- 
ing the truncated length distribution function from the 
family (|7j) and q> \ with Lmax > 20 are indistinguish- 
able from the infinite ranged p{l) which shows that for 
rapidly decaying p(Z)'s the rods with lengths distributed 
in the tail of p{l) don't affect the phase behaviour of the 
polydisperse mixture. In Fig. ^ the values pf, for which 
this transition occur are plotted as a function of the poly- 
dispersity A, for different choices of q using the infinite 
ranged p{l). Also are plotted the cloud isotropic curves 
(the values of p^^' when x — 1) as the solution of the co- 
existence equations (|30| ) and (^l|). For the complete wet- 
ting transition at the wall-isotropic interface the cloud 
isotropic density p^^^ coincide with the wetting density 
value pw As we can see from the figure, both curves are 
practically the same for any distribution function p{l) be- 
longing to the family ( [74| ) with q>l. The general trend 
is the decrease of the transition densities with polydis- 
persity. This is due to the fact that the alignment of 
long rods parallel to the wall is cntropically favored. The 
amount of long rods increases upon increasing polydis- 
persity, therefor this entropic effect increases the total 
first moment parallel to the wall m(0) = mx{Q) + ''nyifi) 
reaching the threshold value for which a second order 
nematic transition appears at the wall. 

Figure ^ also shows the relative stability of the biax- 
ial nematic phase with respect to the wetting transition, 
measured by the magnitude 1 — ph/p^. As we can see, the 
stability does not change too much (from 0.18 to 0.21) 
and in general is an increasing function of A. This means 

u,, is more sen- 




(I) 



that the bulk transition density p 
sitive to the change in polydispersity than the surface 
biaxial transition density p^. 



FIG. 4: In the upper figure the density pi, of the transition to 
the biaxial phase and the wetting density {pw > Pb) are 
plotted as a function of polydispersity A, for q=l, 1.3, 1.5, 1.7, 
2 and 5. The bottom figure shows the relative stability of the 
biaxial phase for the same g's. The points show the maximun 
degree of polydispersity for a given q. In decreasing order of 
A the points correspond to p(/)'s with increasing order of q. 
In all figures the densities, first moments and specie lengths 
are dimensionless (see the text in section ^ . 



B. Bulk and surface phase diagram for slowly 
decaying distribution functions. 

The other selected family of length distribution func- 
tions p(O0(imax — I) with p{l) from ( |76|) and with Ca,q 
and Iq calculated from the usual constraints (the equal- 
ity to one of the zero and first moment). There are ex- 
perimental works with oil-in-water emulsions where the 
droplet size distribution is well described by a log-normal 
distribution function |25[ |. Also was proposed a simple 
statistical model to describe the side distribution of nu- 
clei in films which grow by the absorption of atoms which 
arrive either by impingement or via surface diffusion (ex- 
ample is adsorption of Ag on amorphous carbon [ p6| ) 
The resulting length distribution function is log-normal 
|p7| . In some polydisperse homopolymer systems in a 
single solvent was found untypical topology of the cloud- 
point curves Was concluded on the basis of ex- 
perimental evidence that this behaviour is due to high 
asymmetry in the molecular weight distribution, conclu- 
sion that agree with theoretical calculations using the 
Flory-Huggins theory for polydisperse polymers with log- 
normal distribution function p9[ | . As we can see this type 
of distribution is not rare in polydisperse systems with 
particles which dimensions are determined by some coa- 
lescence process. 

In the set of equations that determine the coexistence 
between the isotropic and nematic polydisperse phases 
there are involved integrals of the type / dll^p{l) exp(/3Z), 
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FIG. 5; Bulk phase diagram of the log-normal polydisperse 
hard rod model with Lmax = 35 (see Eq. |7^ with g = 2). 
The solid lines are the densities of the cloud isotropic (p'^') 
and shadow nematic (p''^') coexisting phases as a function of 
polydispersity (A). The cloud nematic and shadow isotropic 
are the dashed lines. The former exhibit a triple coexistence 
at A = 0.3834 between an isotropic phase and two different 
nematic phases whose densities are marked by black circles in 
the shadow nematic curve. 
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FIG. 6: The first moment distribution function Ip^ {})/po 
for the two coexisting nematics. Continuous lines represent 
the A^i (top) and N2 (bottom) distribution functions while 
the isotropic cloud distribution is plotted with a dashed line. 



with A: = 0, 1. To calculate the first moment or the den- 
sity of the component parallel to the director in the ne- 
matic phase we find integrals of this type with /3 > 0. 
Then for p{l) decaying at large distances slower than 
an exponential these integrals diverge when Lmax ~> oo- 
This means that in the phase diagram the unique sta- 
ble phase is the nematic at any density. When Lmax has 
a finite value there exist a surface phase diagram very 
different for the preceding one, as I will show bellow. 

The bulk phase diagram of the polydisperse hard rod 
system with the length distribution function family j?^ ) 
turns out to be more complex, so I will now describe in 
more details the whole bulk phase diagram and not only 
the cloud isotropic branch (or wetting transition line) , as 
in the preceding case. 

Carrying out the same calculations for coexistence as 
in the previous case, but now using (|7^) with q — 2 
and Lmax = 35, (i.e. the truncated log-normal distri- 
bution function) , we find the bulk phase diagram plotted 
in Fig. ||. In this figure the cloud isotropic-shadow ne- 
matic (a; = 1) and the cloud nematic-shadow isotropic 
(x = 0) densities are shown as a function of A. 

As we can see from the figure there is a polydisper- 
sity value (A = 0.3834) for which in the cloud isotropic- 
shadow nematic curves appear a triple coexistence be- 
tween an isotropic an two different nematic phases (a 
triple point). This behaviour is peculiar of distribution 
functions decaying slower than an exponential for large 
I's and strongly depends on the cutoff Lmax- I have found 
that the two coexisting nematics have different distribu- 
tion functions. The nematic with the highest density 
and first moment {N2) has a bimodal first moment dis- 
tribution function lp^^^{l)/po with the two values of its 
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FIG. 7: The same as Fig. ^, but instead of p the total first 
moment m is plotted as a function of polydispersity. The 
meaning of the different lines is the same as in Fig. U 



maxima having the same order of magnitude, whereas 
the other nematic (Ni) has a unimodal distribution (see 
Fig. I). 

In Fig. 1 have plotted the total moment m of the two 
coexisting isotropic and nematic phases as a function of 
polydispersity. The most interesting feature is the promi- 
nent maximun exhibited by the shadow nematic curve. 
Its existence will be explained below. 

The origin of the second maximum in the function 
lp{l)/po at large polydispersities is very easy to explain 
if we take into account that from Eq. (|3g), the total 
moment of the coexisting nematic phase for a: = 1 is pro- 
portional to the integral of a term which, for large ^'s, 
can be approximated by 
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f{l) = Zexp ( -aln^(-) +4Am? 



(84) 



where Am = m''^' — m^'' > 0. The minimum value of 
this expression is reached at I*, the largest solution to 

nn = 0, i.e. 



2aln(/7;o) 
4Am 



1 



(85) 



If I* < Ly^i^^ we find that the function f{l) has a sec- 
ond maximum at I = imax, apart the other one near 
Iq. From Eq. (|8^) we conclude that increasing polydis- 
persity (decreasing a) I* decreases, so the value /(imax) 
increases as we can see directly from ( ^^ . Also from (|8^ ) 
and the condition I* < Lmax it is clear that f{L„^gj^ is 
an increasing function of Am, thus taking into account 
that m||^^'' < m||^^^ and m||^^^ = rri'j_^\ i.e. the most or- 
dered nematic has less amount of rods perpendicular to 
the nematic director, we find that for 7V2 the second max- 
ima value is higher. For the less ordered nematic (iVi), 
although lp^-^\l)/po can develop a second maximum at 
Lmax, its value is less than the first maximum by several 
orders of magnitude. This is the reason why we call the 
Ni distribution unimodal distribution. 

The maximun in the shadow nematic curve (total mo- 
ment m as a function of A (see Fig. ^)) as function of A 
can be explained again by resorting to Eq. (84). From 
the general cloud isotropic curve behaviour as a function 
of polydispersity Am(> 0) is a decreasing function of 
A. From ( |84| ) the value of f{l) for large Vs depends on 
the two terms in the exponential. For moderate poly- 
dispersities (near the triple point) a rapid decrease in 
a determines the increase of J dlf(l), whereas for large 
polydispersities, where a changes slowly with A and Am 
keeps on decreasing, the second term dominates the be- 
haviour of this integral, which decreases with A. Thus 
at intermediate values m pass trough a maximun. 

When A is larger than its triple point value the stable 
nematic is the bimodal one. The cloud nematic and the 
shadow isotropic curves {x = 0) have the usual behaviour 
with polydispersity, in the sense that they show, for any 
A, a single I-N coexistence, so is expected that for a 
certain x between 1 and this triple coexistence will 
disappear. 

From equation ( p9| ) we can estimate the value x for 
which the contributions of the two coexisting phases 
to the integral have the same order at I = Lmax (the 
length value for the second maxima in the N2 first 
moment distribution function). I obtain x ~ 1 — 
3 exp (— 4AmLmax/3). Taking into account that L^ax = 
35 and that Am is of order of unity we conclude that 
the range of a; at = 1 — a; for which take place the preced- 
ing described behaviour should be of order 10"^". Thus 
the triple coexistence will extend until very small val- 
ues of a; AT. To find the triple equilibrium I have solved 
the set of five equations: Sa = Po't's^, Ta = Po4>rc, 
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FIG. 8: The limited area by the closed curve in the a;jv — A 
plane coincide with the triple coexistence region which finish 
at a tricritical point (the filled square). 



(with a — Ni,N2) and AI1{tn2, SN2) = (see equa- 
tions (^0|), ( pl| ) and (|37|)) whereas po is calculated from 
An(r7Vi , sati ) = 0, i.e. from equation (|3^) with substitu- 
tion N — !■ iVi. The unknowns are Sq, Tq and xjv. The 
definition of E{1) is now 



E{1) = 3(1 - xn) + XN^eiil) + XN2^2il), 



(86) 



where xn = xn-^ + XN2 {xn^ = XN^a, with 9i ~\- 62 = 1) 
is the total fraction of volume occupied by the nematic 
phases and Oi measure the relative fraction of the Ni 
phase with respect to the other. In figure ^ is plotted 
in the xjq — A plane the closed area where the triple 
coexistence was found. Fixing some A in the range 
[0.3834,0.4058] and increasing x^q from zero is reached 
the value for which appear I — Ni — N2 triple coexis- 
tence with infinitesimal amount of the iVi phase {9i = 0). 
Changing the relative fraction of the nematic A'^i with re- 
spect to the N2 we move through the triple coexistence 
line which finish when the nematic N2 fill an infinitesimal 
amount of volume {9i — 1). 

The first moment values for the three phases with xn 
and A inside the triple coexisting region are plotted as 
function of A in figure ^ where the bottom figure is a 
zoom of the top figure for values of m corresponding to 
the coexisting nematics near the tricritical point. The 
shaded area has the same meaning as in figure ||, i.e. 
represent the allowed values of m inside the triple coex- 
isting region. As we can see the gap between the values of 
m corresponding to coexisting nematics decrease mono- 
tonically with A from the triple point corresponding to 
Xn = (the filled circles) until the tricritical point for 
Xn ^ 10~^^ (the filled square). 

In figure arc plotted the first moment distribution 
function lp{l)/po of the three coexisting phases for A = 
0.395 and xn^ — xn- As we can see the phase N2 has 
the second maxima at I < imax. 

I also studied the effect that has the change of the max- 
imun extreme value Lmax on the bulk phase diagram. In 
figure [Tl| I plot the nematic shadow curve for L„iax = 35, 
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FIG. 9: Phase diagram which corresponds to the triple coex- 
istence region of figure^ m^^"' as a function of A (the solid 
line) are the different branches which begin at A = 0.3834 
(the filled circles) and join at the tricritical point (the filled 
square). The dashed line represent the values for mS'K The 
bottom figure is a zoom of the top one around the tricritical 
point. 
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FIG. 11: The shadow nematic curves for different Lmax's. 
In the top: an interval of A near the triple point. In the 
bottom: the whole range of A. In both figures Lmax = 35: 
solid line, Lmax = 30: dashed line, Lmax = 25: dotted line 
and I/max = 20: dot-dashed line. The circles and squares 
represent the coexisting nematic phases at the triple points. 
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FIG. 10: The first moment nematic distribution functions 
lp{l)/po for A = 0.395 corresponding to the triple coexistence 
with = xm- N2: solid line, A'^i: dashed line and /: dotted 
line. 



30, 25 and 20. As we can see the position of the triple 
point moves to larger A's when L,„ax decrease, and also 
there is some Lmax (between 30 and 25) for which dis- 
appear the triple coexistence. Is reasonable to predict 
this behaviour because the second maxima value of the 
shadow nematic {N2) first moment distribution function 
is an increasing function of Lmax and A. Then if L^ax 
decrease the value A needed for keeping the triple coex- 
istence should increase. We conclude also from figure |ll| 
(the bottom figure) that the coexisting densities for A's 
greater than the triple point value decrease with polydis- 
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FIG. 12: The shadow nematic curves for length distribution 
functions with different decaying laws and Lmax = 35. q = 2: 
solid line, q = 2.25: dashed line and q = 2.5: dotted hne. 



persity when Lmax increase. 

The effect of the decaying law of the length distribu- 
tion function for fixed Lmax on the bulk phase diagram 
is similar to the preceding behaviour as we can see from 
figure [1^. I n this case was changed the parameter q in the 
family^Tq) . Increasing q the triple point position moves 
to larger A's and disappear at some value between 2.25 
and 2.5. The second maxima is higher for low decaying 
length distribution functions and then is necessary less 
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FIG. 13: s = mx{0) — my{0) as a function of po for 3 different 
values of polydispersity: A=0.4705, 0.4298 and 0.4112 from 
the left to the right 



amount of polydispersity to exhibit the iVi — N2 transi- 
tion. 

Now I describe the surface phase diagram for the log- 
normal distribution function. Fixing imax — 35 the val- 
ues for the surface parameter s = mx{0) — fny{0) as a 
function of po (the bulk density at infinite distance from 
the hard wall) obtained as the solution of the equations 
( p5| ) and (|4^) for different values of A are plotted in fig- 
ure 



13, For small values of polydispersity (less than those 



shown in the figure) the surface phase behaviour is sim- 
ilar to the characteristic behaviour of the fast decaying 
distribution functions where at some po (in figure this 
value is determined by the intersection of the curve with 
the p axis) appears a biaxial nematic phase through a 
second order transition. Then the order parameter s in- 
crease with density until reaching the complete wetting 
density value (these transition points are marked with 
black circles). 

When A is larger than some value A* w 0.4 the graph- 
ics of s as a function of p develop an hysteresis branch. 
This is an indication of a first order transition between 
two nematic phases at the wall. This behaviour is com- 
mon for all truncated length distribution function with 
slow decay at large Z's and have not analog in the distribu- 
tion function family with fast decay as we have already 
seen. In figure ^ are shown three of these hysteresis 
branches for different A's. The scenarios is the follow- 
ing: when A > A* after the second order transition 
from isotropic in plane phase (ma;(0) = my{0)) to the 



nematic in plane phase (the triangles at figure 13 rep- 
resent the bulk density value for which this transitions 
occur) at some point occur a first order surface transi- 
tion at the wall between two nematics (see the squares 
in figure |l3| ) and then for larger p the complete wetting 
transition (the circles in figure |l^) . The curves in figure 
|l| for A = 0.4112 and 0.4298 represent this behaviour. 
Increasing A the first order nematic- nematic transition 
density value move to the second order isotropic-nematic 
transition value and then for some A'^ these two transi- 



tion points coalesce. For A > A'^ the surface first order 
transition occur between an isotropic and nematic phase 
and then the complete wetting transition as usual (see in 
figure P the curve for A = 0.4705). 

The first order transition points can be calculated ex- 
actly using the local approximation without carrying out 
the full minimization of the grand potential with respect 
to the density profiles. I describe now the followed pro- 
cedure to find them. The excess surface free energy den- 
sity for the non local model which integral over "z" is 
the surface tension (7 = / dzip{z)) can be expressed as 
a function of Ap{z) (Ap(z) = po — piz)), 0{z), r(z) and 
s{z) as 



-po(2r + 9)- 6t. (87) 



= Ap + — 

The gradient of iIj{z) is 



= -Vt(0 + ^po) - \p,p{2z)e-^io'dz'ri.') 



- -PO dlp{l)T * 5; e-^-*'\ 

where I have introduced the notation 

T*6i =t{z- ^)-t(z + ^). 
The gradient of t(z) is related to V9{z) by 



(89) 



(90) 



where 



2 r°° 

VO = -PO dlp{l)e 



2i: + 'dz'r(z') 



- 3P0 dlp{l)e-' ^^^'-''^'^^^'\ (91) 

For the local model (see eqs. ( ^2| ) and (|5^)) the equation 
( ^ ) transforms to 



Vip{z) 



-Pqp{2z) exp (-4t(z)z). 



(92) 



From (|9^) we conclude that the excess energy density 
profile 4'{z) for the local model is a monotonic decreasing 
function of z and that V?/'(0) = 0. In the range of po 
values where there are two solutions (I will call them 1 
and 2) of the set (p5|)-(|4^) (two nematics or one isotropic 
and the other nematic) the more ordered phase has a 
higher value of t(z) near the wall. Then, from ( p^ we 
see that the excess energy density profile of the more 
ordered branch has a smaller slope. Then the transition 
point (the value of po) can be calculated exactly from 
■01(0) = V'2(0) which taking into account the relation 
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FIG. 14: The surface phase diagram of the log-normal poly- 
disperse hard rod fluid with Lmax = 35. Solid line: The com- 
plete wetting line. Dashed line: The first order I-N or N-N at 
the wall. Dotted line: The continuous I-N phase transition. 
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FIG. 15: The zoom of figure |lj. The point is the density 
value of the bulk triple point. The square is the critical point 
of the N-N 2D phase transition. The triangle is the point of 
coalescence of the continuous TN and a first order 1-N 



tp{z) = Ho — n(z) = (where n(z) is the local pressure, 
i.e. the pressure evaluated at p^(z))) becomes in ni(0) = 
112(0). The total integral of ^{z) should be minimized 
to obtain the equilibrium profile, then the length of the 
more ordered phase for higher po 's is calculated from the 
intersection between two branches V'i(^) ^nd '02 (-z)- The 
length of this layer increase continuously from I* as a 
function of po- The value I* depends on the form of 
p{l) for small Vs. For length distribution functions that 
have a continuous decay to zero when / — > 0, P = and 
the transition from less ordered to more ordered density 
profile near the wall is continuous, but if p{l) is zero for 
I < Iq, then I* = Iq and at the transition point coexist 
two profiles with different adsorption coefficients, one of 
them including a layer of length I* of more ordered phase, 
so the transition becomes first order. 

We conclude that although at the transition point 
there is a "2D" first order transition between two differ- 
ent phases, in "3D" the order of this transition depends 
on the form of p{l) for small Vs. The squares in figure 
were calculated using this procedure and the com- 
plete surface phase diagram for the log-normal polydis- 
perse system with Lmax = 35 is shown in figures ^ and 
|l5|. When the polydispersity value is less than the value 
corresponding to the square in figure ^ we find after 
the continuous uniaxial-biaxial nematic transition (the 
transition density is shown with dotted line) the com- 
plete wetting transition (the transition density is shown 
with solid line). When A is between the square and the 
triangle there is also a 2D first order N-N transition in 
between. Finally for A's larger than the position of the 
triangle this 2D first order transition is between isotropic 
and nematic phases. Just at the triangle coalesce a sec- 
ond and first order transitions. The solution of the 
set of equations (^-(^ at z = and the condition 
■01(0) = ■02(0) is equivalent to find the coexistence be- 
tween two phases in the two dimensional bulk system of 
polydisperse hard rods (pz(0) = and rnz(O) = 0) with 



the excess energy density — 2$2i?, where $2i3 i^ 

the excess part of the bulk free energy density of the 2D 
system. The reason why I have found two coexisting ne- 
matics in the bulk phase diagram when xn = only for 
A = 0.3834 (the bulk triple point) whereas for A > 0.4 
exists at the wall a whole range of first order N-N tran- 
sition is the following: The equations (^5|) and (^) for 
z = can be put in the form 



T = |po (^J dllq{l)e-''^ cosh{sl) - 1^ , (93) 
s^'^p^ J dllq{l)e-^'smh{sl), (94) 



which form coincide exactly with the set of equations 
resulting from the equality of chemical potentials be- 
tween the cloud isotropic and shadow nematic phases 
for 2D system with the following variable definitions: 
r = (mi^^)+m(^^))/2-mr\ s = (^r) -™(2^))/2- 
m'^o°\ p[,^^' = 4po/3 and q{l) = p{l)e^P"^'^. Then 
the surface-parent phase from which separate the differ- 
ent coexisting phases have not a log-normal distribution 
form. If the bulk density po is large enough then the 
"surface-parent" length distribution function q[V} has a 
second maxima aX z = Lmax and with a bimodal par- 
ent distribution function is possible to find a nematic- 
nematic coexistence. 

For the non local model is very difficult to gain insight 
about the monotonic behaviour of the density profiles as 
can be seen from the equations (p8|)-(|9l|). In principle 
they can be non monotonic functions of z. Then, the 
location of the transition points between different surface 
phases through a first order transition may differ from the 
local model results. 
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The wall-isotropic and isotropic-nematic 
interfaces 



I have solved the equations (p5|)-(p7|) using the 
Schultz's length distribution function for the bulk 
isotropic phase to obtain the equilibrium wall-isotropic 
profiles. I used the Gauss-Laguerre quadrature to cal- 
culate the integrals over variable I using 10 points and 
using the fast Fourier transform algorithm to calculate 
the convolutions pz * Oi for each I with a grid spacing 
Az = Zniin/10 (?min IS the Smallest quadrature point). 
I have compared the results with those obtained by in- 
creasing the number of points used in the quadrature 
until 20 and the density profiles and the surface tension 
values does not change appreciably for < A < 0.55. 

For polydispersity value A = 0.064 and bulk isotropic 
density pa — 1.02 > ph {pb is the density of the isotropic- 
biaxial nematic transition which is shown as a function 
of A in the inset of figure |l^) I obtain the density pro- 
files Px(z), Py{z) and Pz{z) shown in figure |l6|. In the 
upper figure are plotted the order in plane parameter 
Qb — {Px — Py) / (Px + Py) for different values of po . When 
Pa ^ Pw {pw as a function of A is also plotted in the 
inset of the upper figure) the point z* where separates 
the px and py branches go to infinity continuously as 
should be in the complete wetting transition, with the fi- 
nal profile composed by macroscopic layer of the nematic 
phase. I have checked that jwi = Jwn+Jni for Po = Pw 
(with jwi, IWN and jni the wall-isotropic, wall- nematic 
and nematic-isotropic surface tensions) which should be 
fulfilled at the complete wetting transition. I also have 
found that r*^°) and F^^^ diverges logarithmically with 
coefficients and ai depending on A as we can see from 
figure The general trend is that the prefactor (z/q) of 
the logarithmic divergence of the density adsorption co- 
efficient is a decreasing function (in absolute value) of A, 
whereas the prefactor (ui) of the logarithmic divergence 
of the first moment adsorption coefficient is an increasing 
function (in absolute value) of A. This result is due to 
the fact that the density gap between the wetting density 
Pw — P^^-* and the density of the coexisting nematic phase 
p(^) decrease with polydispersity as we can see from fig- 
ure 18, whereas the moment gap is an increasing function 
of A(see figure ^8|). The surface tension as a function 
of Po is plotted in figure ( p9| ) for three different values of 
A. As I have already shown in section III the behaviour 
of 7 is governed by certain combination of F'-*'-' and F'--'^-' 
(see equation (p^)), thus taking into account that these 
quantities have a logarithmic divergence can be used the 
approximate analytic equation (^^ to estimate the dif- 
ference between the wall-isotropic surface tension "fwi for 
Po < Pw and its wetting value at po = Pw- For compari- 
son I have plotted in figure 19 the exact value of ^wi and 
the approximation ( |37| ) for A = 0.064. If we take into 
the account that the comparison between them should be 
restricted to the narrow neighborhood of pw is surprising 
that they are similar until i = 1 — p/p^, « 0.1, and that 
the maximun position is well estimated. Also is plotted 





FIG. 16: The upper figure: Order parameter profile Qb(z) 
for A = 0.064 and po = 1.05, 1.2, 1.23 and 1.245 (from left 
to right). The inset: Biaxial nematic transition density ph 
and the complete wetting density as a function of A. The 
bottom figure: Density profiles px{z), py{z), pz{z) for po = 
1.245. 
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FIG. 17: The coefficients Vk and ak of the logarithmic diver- 
gence as a function of polydispersity. 



jwi for A — 0.33 and 0.53. 

In figure (the upper one) are shown the isotropic- 
nematic interfaces for po = Pw and A = 0.577. In partic- 
ular are plotted the profiles p{l, z)/ pq{1) for three differ- 
ent Vs. This ratio should tend to 1 at the bulk isotropic 
side of the interface. From figure is clear that the nematic 
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FIG. 18: The quantities 5p 



' (solid line) and 



Sm, = ' — rrv- ' (dashed line) as a function of A. 
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FIG. 19: The wall-isotropic surface tension as a function of 
t — 1~ p/pm for A = 0.064 (solid line), 0.33 (dashed line) and 
0.53 (long line). The dot-dashed line represent the analytic 
approximation (psj) for A = 0.064 



phase has more amount of large rods than the isotropic 
one. This effect is more dramatic for higher A's. In the 
bottom figure I show for fixed / = 1.91 the profiles for 
two different A's. By increasing A the amount of rods 
with this length changes its segregation behaviour. For 
A = 0.333 they are lightly segregated at the isotropic 
phase whereas for A = 0.577 they strongly segregate to 
the nematic phase. 

I have also calculated the wall-isotropic interface in 
the local approximation using the log-normal length dis- 
tribution function for the bulk isotropic phase solving the 
equations (|45|), (^ ) and (|5^). The main reason to use 
the local approximation is that is impossible to select an 
effective quadrature to estimate the integral over / due 
to the specific form of the p{l, z) for large Vs (which can 
develop a second maxima). 

Far from the wetting transition in the region of p^s 
where the slope of ^wi with respect to po is practically 
zero we have calculated the density profiles p{z) and m(z) 
for A = 0.4298. In figure 21 are given these profiles for 
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FIG. 20: The upper figure: the density profile p{l, z) / po{l) of 
the I-N interface for A — 0.577. Are plotted three profiles 
for different Vs. Solid line: / = 1.91, dashed line: I = 1.21 
and dotted line: I — 0.67. The bottom figure: Profiles corre- 
sponding to I = 1.91 for A = 0.577 (solid line) and A = 0.333 
(dashed line). 



po = 0.8397. As we can see from the figure whereas the 
density adsorption coefficient r'"' is negative, the first 
moment adsorption coefficient F^^^ is positive. As I have 
already shown (see eq. (|6^)) the slope of the surface 
tension depends on certain combination of F*^"^ and F*^^^. 
For this particular case the slope is negative although 
F(°) < in contrast with the result obtained using the 
interface Gibbs-Duhem equation for one component sys- 
tem. 

For some A's greater than A* (the critical point value 
for the first order Ni — N2 transition at the wall) were 
calculated the adsorption coefficients F^'^^ and F*^^) as a 
function of po- The results are plotted in figure ^ For 
A = 0.4298 (the solid line in the figure) increasing po 
occur first a second order I — Ni transition and then a 
first order Ni~N2 transition in the plane of the wall. Was 
already shown that the surface length of the N2 phase 
increase from zero resulting in a second order surface 
transition although in 2D is first order. Finally the length 
of the N2 surface phase diverge continuously to infinity 
at the complete wetting transition point. For A — 0.5104 
(the dotted line in the figure) was found a first order I—N 
transition in the plane of the wall and then the complete 
wetting transition. 
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FIG. 21: Density profile p{z) (solid line) and the moment 
profile m{z) (dashed line) for A = 0.4298 and po = 0.8397. 
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FIG. 22: Left figure: The density adsorption coefficient P'"' 
vs. t = l- po/p» for A = 0.5104 (solid line) and A = 0.4298 
(dashed line). Right figure: The first moment adsorption 
coefficient P'^-* vs. t for the same A's. 



As we can see from this figure these surface transitions 
although second order are very hard, i.e. there is not 
a jump in the graphic of r*^°'^^ as a function of t but its 
slopes have a strong discontinuity at the transition points 
which translates in a discontinuity in the second deriva- 
tive of 7vK/ with respect to t as we can see from figure 
p3| . In the figure is also plotted the surface tension for 
A — 0.3837 corresponding to the bulk triple point. For 
this value the transition to the biaxial nematic surface 
phase is second order (see figure [ij) resulting in the con- 
tinuity of (ir^°'^^ /dt (k = 0, 1) at this point and then the 
surface tension curvature is also continuous. 

The coefficients z^o.i and cro.i of the logarithmic diver- 
gence of r*^°-' and r'^' for log-normal p{l) have the same 
qualitative behaviour (see figure |^) as comparing to the 
Schultz p{l), i.e. {vi) is a increasing (decreasing) func- 
tion of A. As I have already shown this behaviour is due 
to the monotonic dependence of the differences p'^^^ — p'-^^ 
(decrease with A) and nipf — mi (increase with delta) cor- 
responding to the coexistence between shadow nematic 
and cloud isotropic phases. As we can see from figure |^ 
m*-''^^ has a strong maximun at A ~ 0.5 so is possible 
that !/i becomes non monotonic function for higher A's. 




FIG. 23: The magnitude A7 = 7 — 7^, vs. t for three difi'er- 





(dotted line). 
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FIG. 24: The coefiicients ft 
vergence of P^'^''^' (see eq 
lated values for A = 0.3837, 
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and (70,1 of the logarithmic di- 
h). The points represent calcu- 
0.4298 and 0.5104. 



D. Enhancement of capillary nematization by 
polydispersity 

Finally I have studied the capillary nematization of the 
slit as a function of polydispersity. In |Q was studied 
this phenomena for the one-component system. Here 1 
pretend to study the role of polydispersity in the possible 
surface transitions governed by confinement. The calcu- 
lations was done in the following way. Fixing the distance 
between two hard walls and fixing the density distribu- 
tion po{l) of the reservoir (which is isotropic) in chemical 
equilibrium with the fluid in the slit, 1 minimize (^9|), this 
time with an external potential which include one more 
wall at z = h. For small values of po the density profile is 
composed approximately by two wall-isotropic interfaces. 
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FIG. 25: The inverse of the width of the slit as a function of 
the reduced density for which occur the capillary nematiza- 
tion. Are plotted three curves for different As: 0.064 (solid 
line), 0.333 (dotted line) and 0.577 (dashed line). 



but there is some (which is less than the complete wet- 
ting density for a single wall-isotropic interface) for 
which the system overcome a first order transition to a 
phase in which the slit is completely filled by the nematic 
phase. This was called capillary nematization in analogy 
with the usual capillary condensation . There is some 
value of the distance between the walls h* (the critical 
point value) for which the transition becomes second or- 
der. The natural question that arise is how depends this 
phenomena with polydispcrsity. Increasing polydispcr- 
sity the entropic ordering effect that has the wall on the 
fluid particles also increase due to the presence of large 
rods. Then is reasonable to predict that the capillary 
nematization will be enhanced by polydispcrsity, i.e. fix- 
ing the distance between the walls the density value po 
for which occur the transition will decrease with polydis- 
pcrsity, and the critical point should move to higher /I's 
as was obtained by numerical minimization. 

The effect that has the polydispersity on capillary 
nematization is shown in figure (p5[), where are plotted 
the transition values h*-t = 1 — pj Pw For fixed t the 
width h increase with polydispersity. 



VI. CONCLUSIONS 

Using simple Zwanzig model in the Onsager limit for 
a polydisperse hard rod system I have study how de- 
pends the phase diagram topology on the selected length 
distribution function p(/). I have found that the relevant 
property which affect the bulk phase behaviour is the de- 
caying law of p{l) at large Vs. For distribution functions 
decaying exponentially or more rapidly was found a single 
coexistence between an isotropic and nematic phase for 
any polydispersity. The cloud isotropic and the shadow 
nematic coexisting densities are decreasing function of A, 
whereas the first moments are increasing functions. For 
distribution functions decaying slower than an exponen- 



tial always we can find a maximun extreme value for the 
lengths of particles (Xmax) for which the bulk phase dia- 
gram turns out to be more complex. Now there is some A 
for which the cloud isotropic and shadow nematic curves 
{xn = 0) develop a triple point: one isotropic (/i) co- 
exist with two different nematics (A^i and N2), one of 
them {N2) with a bimodal distribution function. The re- 
gion of triple coexistence extend until very small values 
of xjv which strongly depend on the decaying law and 
on the maximun extreme value Lmax- For values of a; at 
and A outside the triple coexistence region the shadow 
nematic curve becomes continuous keeping its usual be- 
haviour (a decreasing function of A). Extending in the 
appropriate way the interface Gibbs-Duhem equation for 
the thermodynamic of interfaces to the polydisperse sys- 
tem I have found that the slope of the surface tension 
depends on certain combination of different moment ad- 
sorption coefficients which was confirmed by numerical 
minimization. I have studied also the surface phase di- 
agram of the hard rod system interacting with a hard 
wall for two representative families of distribution func- 
tions. Each of them produce two qualitatively different 
surface phase diagrams. With the particular choice of the 
Schultz's distribution function from the first family were 
calculated the equilibrium profiles and different adsorp- 
tion coefficients. Increasing the bulk density value was 
found some pg for which the system overcomes a second 
order uniaxial-biaxial nematic transition The thickness 
of the new biaxial nematic layer increase from zero and 
diverges continuously at the wetting transition density 
Pw Both Pq and pw are decreasing functions of polydis- 
persity. For the log-normal distribution function (a par- 
ticular choice from the second family) the surface phase 
diagram changes dramatically. The second order I — Ni 
transition at the wall is followed by a two dimensional 
first order transition between two nematics (A^i and A^2). 
The thickness of the N2 layer increase from zero contin- 
uously and diverge at the wetting density. Increasing A 
this A^i — N2 transition coalesce with the second order 
/ — A^i transition and for higher values we find a single 
first order I—N2 at the wall followed by the complete wet- 
ting. For the log-normal distribution function I have ob- 
tained using the local approximation that although very 
hard (the slopes of the adsorption coefficients at the tran- 
sition points have strong discontinuities) these transitions 
are second order in nature. Was studied also the segre- 
gation phenomena occurring at the isotropic-nematic in- 
terfaces. The numerical calculations show a preferential 
segregation of long rods to the nematic side of the inter- 
face being very sensitive this phenomena to the degree 
of polydispersity. Finally was shown that polydispersity 
enhance the capillary nematization of the slit. 
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